Primordial helium recombination III: Thomson scattering, isotope shifts, and 

cumulative results 
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Upcoming precision measurements of the temperature anisotropy of the cosmic microwave back- 
ground (CMB) at high multipoles will need to be complemented by a more complete understanding 
of recombination, which determines the damping of anisotropies on these scales. This is the third 
in a series of papers describing an accurate theory of He I and He ll recombination. Here we 
describe the effect of Thomson scattering, the ^He isotope shift, the contribution of rare decays, 
colhsional processes, and peculiar motion. These effects are found to be negligible: Thomson and 
''He scattering modify the free electron fraction Xe at the level of several x 10 The uncertainty 
in the 2'^P° — 1^5* rate is significant, and for conservative estimates gives uncertainties in Xe of 
order 10~^. We describe several convergence tests for the atomic level code and its inputs, derive 
an overall Ce error budget, and relate shifts in Xe{z) to the changes in Ce, which are at the level of 
0.5% at £ = 3000. Finally, we summarize the main corrections developed thus far. The remaining 
uncertainty from known effects is ~ 0.3% in Xs- 

PACS numbers: 98.70.Vc, 95.30.Jx 



I. INTRODUCTION 

Cosmological recombination determines the evolution 
of the free electron fraction as the universe becomes cool 
enough for bound atoms to form. Concurrent with this 
is a (cosmologically) rapid drop in the Thomson opacity 
that decouples the motion of the photons from that of the 
baryons; the perturbations in the photons are observable 
today as the anisotropies of the cosmic microwave back- 
ground (CMB). Thus, many properties of the measured 
CMB anisotropy are strongly dependent on the free elec- 
tron history. The free electron fraction sets the damping 
scale for temperature anisotropies in CMB through Silk 
damping, and contributes to the acoustic damping scale 
P, U, S @1 ■ These are manifest in the positions of 
the acoustic peaks in the temperature anisotropy power 
spectrum, and by the suppression of the anisotropy power 
on small scales. Because of the large number of modes 
on the sky, the temperature anisotropy power can be 
measured very accurately, even with small survey cov- 
erage. These modes are the target of a new generation 
of high-precision, small-scale temperature anisotropy ex- 
periments 0, i, i, M M m, M M, M, MM,m, and 

will provide stronger constraints on baryonic and mat- 
ter fractions, the primordial spectral slope ris, and its 
possible scale dependence. 

The physics underlying recombination is well- 
developed and the physical circumstances are simple: the 
universe is homogeneous, there are no elements heavier 



* Electronic address: 
t Electronic address: 



switzer@princeton.edu| 



chirata@sns.ias.edul 



than Li, and all material is in the gas phase. The sub- 
tleties in the recombination history emerge from both 
the very high accuracy required, and the large num- 
ber of processes contributing. Here, rare processes are 
important because many of the fast allowed processes 
have reverse processes that quickly come to equilibrium, 
blocking overall progress of the forward direction. Exam- 
ples are recombinations directly to the ground state and 
2^P° l^S decay in He i: these are fast but access only 
certain regions of photon phase space {hi/ > 24.6 eV and 
the 21.2 eV line, respectively). Instead of bringing the 
ionization state of hydrogen or helium into Saha equilib- 
rium, these processes merely boost the number of pho- 
tons in certain regions of phase space to values far greater 
than the Planck distribution Af = i/(e''^/'^BT _ -^^ would 
predict - until some process (such as Sobolev escape in 
the case of the 21.2 eV line) removes the photons. It 
is difficult to have a prior notion of what rare processes 
will be significant, and to what extent. In general, pro- 
cesses that modify the overall recombination history ei- 
ther break the equilibrium in allowed lines by removing 
photons, or provide an entirely independent path to the 
ground state. Examples are H I photoionization opacity 
and two-photon decays, respectively. 

Recently, several corrections to the theory of He I re- 
combination have been proposed [l^, [13, [111, [H]. This 
is the third of a series of papers describing new effects 
and remaining uncertainties in He i recombination. Pa- 
per I (Switzer and Hirata astro-ph/0702143 ) describes 
the base recombination model, the influence of the feed- 
back of spectral distortions between lines, and the effect 
of continuous opacity from H i photoionization in trans- 
port phenomena during He i recombination. Paper II 
(Hirata and Switzer ,astro-ph / 0702 144| ) describes absorp- 



2 



tion of non-thermal radiation in two-photon processes 
from n = 2, non-resonant two-photon decay processes 
from n > 2, and the effect of finite resonance linewidth. 
(Henceforth we refer to these as Paper I and Paper II.) 
Here, we describe several additional effects that are negli- 
gible to recombination, summarize the overall magnitude 
and convergence of effects studied, and give an error bud- 
get in free electron fraction and the anisotropy power in 
the CMB. 

Because of the large thermal velocities of electrons 
in the primordial plasma, Thomson scattering results 
in a very wide redistribution of energies. A typical 
Thomson scattering will shift the photon's frequency by 
~ 4 THz. For comparison, the Dopplcr width associ- 
ated with scattering through atomic transitions in He i 
is only ~ 50 GHz. The redistribution width has impor- 
tant consequences for radiative transport within the line 
subject to Thomson scattering. We will see that in the 
case of the intercombination lines {n'^P° — 1^5* in He i), 
it is unlikely for a photon to scatter with an electron as 
it traverses the width of the line that is optically thick to 
incoherent processes because of the low differential opti- 
cal depth [dTe/dv) to Thomson scattering. However, a 
photon can be scattered by an electron as it redshifts far 
from line center - as it redshifts, the trajectory integrates 
sufficient differential depth for Thomson scattering to be 
likely. Because of the large typical frequency shift from 
electron scattering, this exchange can re-inject the pho- 
ton onto the blue side of the line, where it has a high 
probability of redshifting back into the line and being 
absorbed. This will further excite He i atoms and inhibit 
recombination overall. The allowed lines have sufficient 
natural width for Thomson scattering to pull photons 
from the region of the line that is optically thick to in- 
coherent processes, which will tend to relax it and accel- 
eration He I recombination. Thus in 2^P° — l^S", which 
is most significant for He i recombination, the effect of 
Thomson scattering in the far red wings opposes the ef- 
fect of scattering closer to line center. 

Thomson scattering is a subtle effect because it is tied 
to the radiation profile on the red side of the line, which is 
determined by scattering in the wings and H i continuous 
opacity. It is addressed here by including it in a photon 
Monte Carlo from Paper I. We will sec, further, that the 
effect of Thomson scattering is highly suppressed when 
we include feedback between lines. In the final analysis, 
including it docs not significantly alter the recombination 
history. 

The second main effect examined here is the isotope 
shift between '^He and '*Hc, dominated by the nuclear 
mass difference. The isotope shift between '^He and '*Hc 
2^P°-l^S is j/(4He)-j/(3He)=263 GHz. Photons on the 
red side of *He i 2^P° — \^S can scatter off of a "^He atom 
and be moved farther into the red damping wing where 
they are more likely to escape. However, in the ordinary 
Sobolev theory, this effect is cancelled (on average) by 
the reverse process where '^Hc scatters a photon blue- 
ward. This is because the radiation phase space density 



on the red side of an optically thick line in an expanding 
background is flat (A/"=constant) to a very good approx- 
imation [2^, [23, [23 • As we saw in Paper I, continuous 
opacity from H i produces a gradient in the radiation 
phase space density, as progressively more photons that 
escape the line arc absorbed, moving redward. Thus, 
once the H i population becomes significant at z < 2200, 
the radiation phase space density on the red side of the 
line is no longer flat. This means that, on average, more 
photons can be scattered redward (because the "pool" 
is larger) than blueward, and the escape probability in- 
creases. We find that this causes a small modification to 
the free electron fraction, |Axe| < 2.5 x 10~^, peaking at 
z - 1900. 

In Paper I we calculate the modification of rates in the 
■n?P° — and ri^D — l^S series by continuous opac- 
ity. Here we describe additional sources of error in the 
recombination calculation related to rare processes and a 
breakdown of the relative contributions of the rare pro- 
cesses. There is considerable disagreement on the spon- 
taneous 2'^P° - l^S rate in the literature [H, [13, TH]. 
This is a significant systematic overall in He i recombi- 
nation, resulting in a maximum uncertainty at 2; ~ 1900, 
of roughly |Axe| < 1 x 10~^. The n^D — series and 
n^po — 1^5' for 71 > 3 are found to give negligible modi- 
fications to the He i recombination history. 

In contrast to the situation in the interstellar medium, 
collisional processes during cosmological recombination 
are subdominant because of the low baryon/photon ra- 
tio. During He i recombination, the highly excited states 
are also so close to equilibrium that electron collision- 
induced transitions between bound and free states cannot 
do anything (Sec. ITVB|) . Two remaining possibilities for 
collisional processes that modify the He I recombination 
history are charge transfer reactions with H i (through a 
collision, a neutral hydrogen atom transfers its electron to 
singly-ionized helium atom) and collisional de-excitations 
to the ground state (where collisions de-excite and atom 
to its ground state without the emission of photon that 
will further excite that transition) [2^, [s^l • In general, 
charge transfer requires an H i population. Continuous 
opacity (described in Paper I) causes He i recombination 
to finish as soon as even a small amount of H i is present, 
and we will see that collisional charge transfer plays a 
subdominant role. In Sec. IIVBI we also show that the 
collisional de-excitations to the ground state are negligi- 
ble from the combination of temperature, free electron 
fraction and collision strengths. 

Throughout the series of papers we consider transport 
processes during He i recombination that depend explic- 
itly on the radiation profile around the line and the line 
profile. This is in contrast to Sobolev methods for the 
escape probability, which do not depend on the line pro- 
file. A significant concern is that we have ignored the 
correlated, peculiar motions in the gas. In astrophysi- 
cal systems, microturbulance leads to broadening of the 
line profile in addition to the broadening from thermal 
motion. Indeed, correlated motions in the recombination 
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plasma arc of the same order as the thermal velocities. 
In Sec. IIV CI we argue that regions of correlated motion 
arc much larger than the distances a photon traverses 
as it scatters through the line. Thus, to a very good 
approximation, the recombination picture is unchanged. 

We propose three main corrections to the standard He i 
recombination model based on the studies in this series 
of papers: (1) inclusion of feedback between allowed lines 
from redshifting of the spectral distortion from n + 1 to 
n, (2) inclusion of rates from intercombination processes 
[20|, and (3) treatment of H i continuous opacity in trans- 
port within 2^P°-1^S and 2^P°-1^S. These corrections 
were addressed in Paper I, and the results from Paper II 
and this paper have little effect. (Note that, throughout, 
we have performed very detailed analyses of the effects 
that we find to be important, or in cases where only a 
detailed analysis can determine whether the effect is im- 
portant. Conversely, we have not investigated processes 
in detail if a simple ordcr-of-magnitude argument shows 
that they are unimportant. The approximate treatments 
are indicated with a ~ sign in front of the magnitude of 
the effect in the summary in Table |l| . 

In this paper we break the discussion of effects into 
several sections. In Sec. |lll we describe the effect of 
Thomson scattering in the n = 2 allowed and intercom- 
bination lines. In Sec. IIIII we treat the isotope shift 
of ^He relative to ■^Hc in the transition 2^P° - l^S. 
Next, in Sec. llVi we study several additional uncer- 
tainties: treatment of intercombination and quadrupole 
lines (Sec. IIV A[) : collisional processes (Sec. IIV Bp : pe- 
culiar velocities (Sec. IIV CP : and convergence of the nu- 
merical code (Sec. IIV Ep . In Sec. |Vl we relate modi- 
fications in the He i recombination history to modifica- 
tions in the temperature-temperature (TT), polarization- 
temperature (TE), and polarization-polarization (EE) 
anisotropy power as a function of multipole. In Sec. IVII 
we summarize the contributions to Xe{z) and anisotropy 
error budget associated with effects investigated in this 
series and describe the prospects for including these ef- 
fects in a fast multilevel recombination code for CMB 
codes. Appendix [X] describes an explicit form for the 
electron scattering kernel, and Appendix|B]describes cos- 
mological perturbation theory relevant for the helium pe- 
culiar velocity. 



II. THOMSON SCATTERING 

Consideration of the effect of Thomson scattering 
in line transport during He i recombination is well- 
motivated because the bulk of electrons have not rccom- 
bined, and their large thermal velocity dispersion relative 
to "^He atoms yields new transport behavior. In this sec- 
tion we analyze this behavior and relevance to the He i 
recombination history. Electron scattering should also be 
folded into radiative transport in H i recombination; this 
is beyond the scope of discussion here and is deferred to 
later work. This is a harder problem because of the high 



optical depth, broad linewidth and partial redistribution 
in the H i 2p <-> Is system, as well as the more stringent 
accuracy requirements. In particular the very high op- 
tical depth (of order 10^) may require a Fokker-Planck 
or hybrid approach to the problem instead of the Monte 
Carlo methods applied here. 

As in Paper I, a photon Monte Carlo will be the 
workhorse for calculating the modified escape probability 
in the He i n^P° - 1^5, n^D- \^S, and 2^P° - 1^5 lines. 
This is described in Sec. IIIBI The physics of the scat- 
tering kernel is described in Sec. Ill CI and Appendix [A] 
We also develop analytic methods (Sec. IIIPp based on 
this kernel to check the Monte Carlo with forbidden lines 
in restricted cases. In forbidden lines, the situation is 
simplified relative to 2^P°— because the photons are 
completely distributed across a line which is narrow com- 
pared to the thermal width of electrons, so that accurate 
analytic transport solutions are feasible. We find that 
Thomson scattering is negligible overall due to cancel- 
lation of several competing effects and suppression from 
feedback of spectral distortion between lines. 

A. Physical setting 

In Paper I, we considered line radiative transport in a 
homogeneous gas on an expanding background subject to 
coherent scattering through the line, incoherent processes 
in the line, and H I photoionization. Here, we extend this 
picture to include Thomson scattering from a thermal 
electron distribution with temperature T,„ (though this 
is nearly identical to during He i recombination) . De- 
pending on the electron's velocity, a photon will be widely 
redistributed because the characteristic scattering width 
for a photon of frequency vq is ~ z/q \/k^T^-al (nieC^)- This 
redistribution is shown at z = 2500 in Fig.[Tl and is much 
larger than the Dopplcr width of the ^He line, ~ 50 GHz. 
Fundamentally, the only reason for solving the transport 
problem is to find the probability Pose that a photon 
emitted by an atom will be absorbed through incoherent 
processes in another atom before it can escape the line. 
That is, the full matrix of probabilities (that a photon 
which has just been emitted or scattered will be emitted 
or scattered by given process before it escapes) is imma- 
terial, aside from its contribution to Pose, as in Paper 
I. We will consider the modification to escape probabil- 
ities from Thomson scattering in n^P° — 1^5* (n < 6), 
n^P" - l^S (n < 5), and n^D - l^S (n < 6). 

Electron scattering shares some, but not all of the fea- 
tures of continuous opacity from neutral hydrogen, devel- 
oped in Paper I. Like continuous opacity in H i electron 
scattering presents a differential optical depth approxi- 
mately flat in frequency, 

Electron scattering would naively be expected to become 
important if ry~^ starts to fall within frequency width of 
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FIG. 1: The electron scattering kernel at z = 2500 spans sev- 
eral THz, allowing photons in the far red side of the line to 
be scattered to frequencies significantly above the line center 
frequency. Rather than escaping on the red side of the line, as 
in ordinary Sobolev escape, the photon can then be absorbed 
(and likely be re-emitted) by incoherent processes in the line 
before it can escape, thus reducing its escape probability. Be- 
cause the radiation phase space density is higher on the red 
side of the line, a photon is more likely to scatter from below 
to above the line frequency, thus decreasing the overall escape 
probability. Once continuum opacity becomes significant in 
the allowed lines, this can also remove trapped photons. 



the line that is optically thick to incoherent processes. 
Airline- In Paper I, this allowed a small H i population 
to remove photons from the line core, so that they can 
no longer further excite the line; thus increasing the es- 
cape probability. The relevant comparison for electron 
scattering is plotted in Fig. [21 which makes it clear that 
electron scattering acts only on frequency scales an order 
of magnitude larger than the incoherent optically thick 
width of the line. Therefore one might expect an accel- 
erating effect of order ~ 10% due to electron scattering. 
In reality this number will be smaller than ~ 10% be- 
cause the Thomson opacity - unlike H i photoionization 
- is from scattering, not absorption. Some of the scat- 
tered photons will re-enter the line, and some photons 
that are on the red side of the line and would otherwise 
have escaped will be scattered into the line, or to its blue 
side (where they redshift back into the line). If there are 
more photons on the red side of the line than the blue, 
more photons are scattered bluewards than redwards on 
average, and the escape probability can be decreased. In 
particular, this means that it is also possible for Thomson 
scattering to delay He i recombination. 



B. Monte Carlo simulation of line transport with 
Thomson scattering 

In this section, we use a Monte Carlo method to es- 
timate the escape probability for a photon trapped in 
He I lines. For reasons described in Paper I, a Monte 
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FIG. 2; The inverse differential optical depth due to Thomson 
scattering of free electrons {r;^^), compared to the inverse dif- 
ferential depth to H I photoionization (??c^), the line Doppler 
width {Auo), and the optically thick linewidth to incoherent 
processes through 2^P°-1^S. Here we can see that a photon 
trapped in a line is relatively unlikely to interact with an elec- 
tron before it escapes. For z < 2500, electron scattering acts 
on scales only roughly an order of magnitude larger than the 
optically thick linewidth to incoherent processes. The behav- 
ior has three regimes for 2^P°-1^S: (1) for z > 2200 electron 
scattering infiuences transport by acting over large scales, re- 
tarding recombination by injecting photons to optically thick 
regions on average; (2) H I continuous opacity suppresses the 
phase space density on large scales: when electron scatter- 
ing occurs, it is much closer to the line's center and removes 
trapped photons, accelerating He I recombination; and (3) 
z < 1900 where H I continuous opacity acts within the line 
and dominates most of the behavior. 

Carlo is a useful procedure to account for the numer- 
ous effects active in the radiative transport within lines, 
despite the computational cost. In particular, the large 
width of the electron scattering kernel compared to the 
width of the He i lines suggests that a Fokker-Planck 
approach is not well suited to evaluating electron scat- 
tering. (The Fokker-Planck operator treats transport of 
numerous scattering events as a continuous process, and 
is accurate if the line profile varies slowly over this range. 
In electron scattering, a photon in the far red can be 
scattered to blue, and this assumption is invalid. We will 
consider an alternative analytical approach in Sec. Ill Dl 
that is useful for testing the Monte Carlo.) 

Kinematically, the frequency shift in a photon scatter- 
ing from an electron in the non-relativistic limit is 

Az/ = (/|[ + a)(l - cosx) - /j.sinx, (2) 

where f = VQwjc is thermally distributed (|| and _L de- 
note the components parallel to the direction of propaga- 
tion of the incident photon, and perpendicular to it but 
in the plane of scattering, respectively), for the photon 
frequency vq, a = hv^ j (me,(?\ and x is the scattering 
angle between the direction of propagation of the incom- 
ing/outgoing photons. The recoil term a is from differ- 
ences in the electron kinetic energy and can be dropped 
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(it is of order 200 GHz, small compared to a frequency 
shift of several THz typical in the exchange due to the 
Doppler shift). This is nearly identical to the case of 
coherent line scattering, except that the scattering cross 
section is flat in frequency and struck atom's velocity is 
not conditioned on the incoming photon frequency. This 
can be trivially added to the Monte Carlo developed in 
Paper I by adding an additional scattering process with 
differential depth r/e, from Eq. ((T|). 

As in Paper I, we simulate the escape probability over a 
grid in {a;Hei, z} with continuous opacity derived from H i 
populations in Saha equilibrium (accurate until the end 
of He I recombination, where it ceases to matter anyway, 
once He i has almost fully rccombincd) . The grid is taken 
over 11 linearly-spaced points in rcdshift spanning z ~ 
1400 to z = 3000, and 21-logarithmically spaced points 
in XHci from 2 x 10~^ to 0.08. This is log-interpolated 
along xhcI and Pcsc, and linearly interpolated along z in 
the level code. 

The results of the Monte Carlo calculation of the 
escape probability for a typical recombination history 
are shown in Fig. [J] for 2^P° — l^S and in Fig. [5] for 
2^P° — l^S. (Fig. d] also has results from an analytic 
method developed in Sec. IIIDI for the escape probabili- 
ties in 23p° - l^S for z > 2000. Fig. [5] includes Monte 
Carlo results for ■^He scattering, developed in Sec. Illlh 
Fig. [5] indicates that Thomson scattering leads to a re- 
duced escape probability for z > 2100 in 2'^P° — l^S" and 
2^P° — l^S, as photons that redshift out of the line are 
Thomson scatttered back on onto it, where they are less 
likely to escape. Once H i opacity becomes important, 
fewer photons are able to travel far enough to be Thom- 
son scattered, and are, instead, more likely to photoionizc 
an H I atom. 



C. The electron scattering kernel 



The electron scattering redistribution kernel has accu- 
rate approximate expressions over a wide range of phys- 
ical scales In the non-rclativistic limit considered 
here, the Thomson scattering kernel with dipolar angu- 
lar distribution (s^ is adequate. We will actually use 
both the kernel and its characteristic function, which 
can be obtained in closed-form from the angular depen- 
dence (Eq. [2) as follows. The electron velocity compo- 
nents /|| and /_L are normally distributed with variance 
'^D = '^oknTni/i'mec'^). Thus, for a fixed scattering angle 
X, the change in frequency is equal to the sum of two 
normal distributions centered around zero, and we write 
the variance of the photon frequency for a given 



Integration against the angular redistribution probability 
gives the redistribution kernel 



P{Av) = ( -j= — — exp 



(4) 



The dipole angular redistribution for electron scattering 
is 

3 

-Pdipoio (x) f^X = g (1 + cos^ X) sin x dx- (5) 

The kernel exists in the literature [s^ in terms of Gaus- 
sian and error functions. In Fig. [1] we compare the kernel 
(in physical frequencies Av) for dipole scattering with a 
with a Monte Carlo calculation. For the analytic work 
of Sec. Ill Dl (needed to test the Monte Carlo) it is more 
convenient to work in the Fourier domain, so instead of 
P{Aiy) [s^l we will concentrate on the characteristic func- 
tion of the photon Ai> distribution. 



vj{k) 



ikA 



^)ai' — 



ikA 



"PiAi') dAv. (6) 



This has a closed form solution for the dipole angular 
redistribution function that is given Appendix [Xj 



D. An approximation to transport in 
Doppler-width dominated lines with electron 
scattering 

In the limit that the linewidth is small compared to 
the characteristic redistribution width for electron scat- 
tering and the H i photoionization opacity within the line 
(r/cAi/iino) is small, we can derive an approximate solu- 
tion for the modification to the escape probability with 
complete redistribution. This is a reasonable approxima- 
tion in the case of the intercombination lines in He i re- 
combination, where the characteristic width is k, 10^ GHz 
and the electron scattering redistribution width is several 
THz. 

To find the escape probability using an analytic 
method, one generally solves for the difference between 
the phase space density of radiation in equilibrium with 
the line and the actual radiation phase space density in- 
tegrated across the line profile, A/l — M (see Paper I). 
Taking Au relative to the line center (y = Av -I- t'linc), 
the transport equation with electron scattering with 
complete redistribution, continuous opacity and electron 
scattering is 



d{Av) 



meC^ 



- cosx). 



(3) 



r;c[AA(A:/) - Ab] + rs0(t.)[AA(Aiy) - AL 



N{Aiy) 



- j p{Av- Au')U{Av')dAv' . {7) 
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The approximation we take here is that (j){l^v) (5(A;/). 
which reduces this equation to 



where 



dN 



+ V. 



- j p{Av - Av')N{Av')dAu' 



where 



I = N{-e)-N{+e) 



(8) 



(9) 



is the jump across the Une. If we consider Eq. ([7]) in 
the immediate vicinity of the hne, we can see that the 
right-hand side contains a delta function at Av = 0, and 
hence there is a jump in M at this value. The jump con- 
dition is obtained by considering the integral <f>(A;^) ~ 



(t){Av')dl\v', which varies from at Av = 
at Av = +e. The relevant terms in Eq. ([5]) are 



dJ\f 



so that Af — A/l oc e"^^*. This has the solution 



-e to 1 



(10) 



(11) 



and the phase-space density averaged across the line is 

^ = f {-^L - Wl - AA(+e)]e"^^(i-*) I d$ = - 

(12) 

In order to proceed further, in particular to evaluate /, 
we need one more relation between JV{—e) and M{+e). 
This can be obtained by transforming Eq. ([8]) into the 
Fourier domain: 

ikC{k) = ricC{k) + i]e [1 - w{k)]C{k) - /, (13) 

where 

C{k) = J [U{Aiy) - Nc] e*'=(^'')d(Ai/). (14) 
This has the simple solution 

"'^' = „c + ,Mi-^Wl-.^ - 

The phase space density averaged over the jump is then 
the inverse Fourier transform, 



C{k)dk, (16) 



i ZTT 



where PP denotes the principal part. Using Eq. (jlSp . 
this can be re-written as 



-PP 



77c + ?7e[l - n7(/c)] 



-dk 



ttJo k^ + {ric+Ve[l-w{kW'^''- ^^^^ 

We can solve for / by algebraically combining Eq. pT]) 
with the definition of / (Eq.^ and with Eq. (HI]). This 
gives 



2 



q + coth(rs/2) ' 



l + q + 

= (Ml-Mc 

From Eq. ((T2]) we find 

AAl-ATs. -(AAl-AAc) ^ I. 

giving the modified escape probability 

2 



(19) 



(20) 



(21) 



AA(+e)+AA(-e) = 2AAc+/g, 



(17) 



"^^^ rs[g-fcoth(Ts/2)]- 

Taking the limit that r/e ^ 0, one can sec that q —> 1, re- 
gardless of the continuum opacity, rjc . Through the use 
of hyperbolic function identities the q ~ I case can be 
shown to give the Sobolev result, Ppsc = (1 ^ e^'^^)/'''s- 
This is by construction of the approximation - for an in- 
definitely thin line, finite continuous opacity can not af- 
fect transport "within" the line. We evaluate Eq. (fT8|) for 
q using a 15-point Gauss-Kronrod rule, using the closed 
form solution for zu^k) from Appendix [Xl The value q is 
shown as a function of rye for several choices of the con- 
tinuum differential optical depth in Fig. [3l We can also 
use Eq. (PT|) to find the modified escape probability for a 
typical recombination history and compare with Monte 
Carlo results. This is shown in Fig. H) 

Note that the analytic method derived here is only 
applicable when continuous opacity does not act within 
the optically thick width of the line. In principle, if the 
intercombination lines remain narrow compared to r]^^ 
until z < 1800, then He i will have nearly finished re- 
combination and the analytic method of Eq. (pij) would 
be appropriate (even though it is technically incorrect af- 
ter that point, it ceases to matter). In a comparison to 
Monte Carlo results in Fig. [H we see that there is signif- 
icant departure from the assumptions built into Eq. (j2ip 
starting at z '--^ 2000. Thus in the level code we use an 
interpolated grid of probabilities from the Monte Carlo, 
and the method here should only lend confidence to the 
Monte Carlo result for z > 2000, where electron scatter- 
ing matters most. 

III. ^He AND THE ISOTOPE SHIFT 

Up to now, wc have considered the He i lines to have 
Voigt profiles. This is not quite correct, because the tran- 
sition frequencies in '^He are shifted slightly to the red of 
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FIG. 3: The quantity q—1 (defined in Eq. (fT8|l ). which quan- 
tifies the departure from the Sobolev theory due to electron 
scattering as a function of the electron scattering differen- 
tial opacity, for several values of the continuum depth for 



2^P°-liS at 



2500, which sets the electron temperature. 
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FIG. 4: A comparison of analytic and Monte Carlo treat- 
ments of transport in 2^P° — l^S which give a modified es- 
cape probability due to electron scattering. This particular 
trajectory of xuci is for the recombination history developed 
in Paper I. The analytic description of electron scattering in 
narrow lines with complete redistribution of Eq. (|21|l agrees 
well with Monte Carlo results at early times. Once continuous 
opacity becomes important within the line, the Monte Carlo 
gives a dramatic increase in the escape probability. Built in 
to the analytic method, however is the assumption that the 
line is indefinitely narrow. Thus, the analytic method breaks 
down z < 2100 and approaches the Sobolev theory. (The 
history without electron scattering is also consistent with the 
analytic method for complete redistribution and continuous 
opacity developed in Paper I.) 



than analogous transition in "^Hc. In the usual Sobolev 
approximation this does not have any effect because the 
hue profile is irrelevant. However, when H i opacity is in- 
cluded, it is possible that splitting the line into two pieces 
(one for '^Hc and one for ^Hc) would speed up recombi- 
nation. This could occur cither if photons rcdshifting out 



FIG. 5: The fractional modification to the 2^P° - l^S escape 
probability due to electron scattering and ^He scattering. The 
effect of electron scattering can be split into two regimes: 1) 
for z > 2300 where most of the effect from electron scattering 
is due to modification to transport far in the wings and the 
escape probability is decreased, 2) for z < 2300, where much 
of radiation in the red wing is absorbed through H I photoion- 
ization, and electron scattering begins to eject more photons 
from within the line, on average. Thus, for 1600 < z < 2300 
electron scattering slightly increases the escape probability. 
^He can be seen to be a small effect-also note that both ef- 
fects are at the level of the noise in the MC measurement. 
These required 3 days across 50 x 3 GHz nodes, so improved 
statistics would require significant computing time. 



of the ''^He line are absorbed by H I before reaching the 
■^Hc line, or if scattering in the '^He line can transport 
photons farther into the red wing of the ''He line where 
they are more likely to escape. In practice the effect of 
■^Hc in radiative transport is more subtle than these sim- 
ple arguments because the '^He line overlaps both the red 
damping tail and the Dopplcr core of "^He. This section 
describes how '^Hc scattering is included in the photon 
Monte Carlo simulation. We argue that the overall con- 
tribution to the escape probability is negligible. 

Among the n^P° - l^S series, the 2^P° - l^S line is 
the primary contribution to the He i recombination rate. 
We will consider the modification of radiative transport 
in 2^P° — l^S caused by ■^He scattering photons from 
the main ^Hc line. If the modification to 2^P° — l^S is 
negligible, then it may be assumed that the higher or- 
der n^P° — l^S series lines have lesser contribution to 
He I recombination. (Even if "^He led to a larger mod- 
ification in the escape probability for n > 2, the levels 
arc so sparsely populated that their overall contribution 
is subdominant.) The isotope shift is inconsequential to 
transport in intercombination and quadrupole lines be- 
cause the analogous '^Hc feature is very optically thin. 

The 2ip° - l^S" line has an isotope shift [H of 



i/(^He) - i/C^He) = -263 GHz. 



(22) 



This difference is approximately four thermal widths of 
the "^He line during most of He i recombination, so both 
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the separation of the hnes and the larger Doppler width 
associated with '^He must be considered. Technically the 
2^P° level in "^He also has hyperfinc structure, with a 
sphtting of 20.8 MHz [H between F = 1/2 and !f = 3/2 
levels, but this is well inside the Doppler or natural width 
of the ■^He line and so can be neglected. Throughout, 
we use an abundance ratio by number of ■^He:^He of 
/4(3He) = 1.1 X IO-VO.079 = 1.4 x lO""^ [H. 

Reciprocity holds to a very good approximation in "^He 
scattering, i.e. the rate for scattering a photon from fre- 
quency vi to 1^2 is equal to the reverse rate from 1^2 to vi . 
(As discussed in Paper II, this is violated if hi/i — hi^2 is of 
order A:BTm, but for the linewidths relevant in helium re- 
combination this does not happen.) For a radiation field 
whose phase space density is constant across the width 
of the "^He line, scattering by '^He cannot have any ef- 
fect. On average, for each photon scattered to the red 
side of the '^He line, there is one scattered to the blue 
side of the '^He line. In He i recombination without con- 
tinuous opacity from H i, the radiation on the red side 
of the *He 2^P° — l^S line is very flat over the width 
of the ^He hne [H, [13, [11. Similarly if the ^He and 
^He level occupation probabilities are equal and there 
is no H I opacity, then the radiation phase space den- 
sity produced in the optically thick ^He line is already in 
equilibrium with the excitation temperature of the "^Hc. 
Therefore without H i continuum opacity the '^He can 
have no effect. However, continuum opacity from H i be- 
comes significant for z < 2200, and produces a gradient 
in the radiation phase space density on the red side of 
^He 2^P° — 1^5. The interaction of partial redistribution 
through ^He and '^He (though only coherent scattering 
in 3 He is significant) and H i continuous opacity is best 
dealt with numerically. 

The Monte Carlo provides a convenient way to include 
scattering through ^He in the radiative transport in the 
^He 2^P°— 1^5* line. Because we take an identical occupa- 
tion history for ^He and *He, and the isotope shift is small 
compared to line-center frequency, the optical depth is 
well-approximated as T5(3He) ~ /4(^He)T5('''He). The 
recoil kinetics and the Doppler width are both modified 
by the mass difference between '^He and ''He. To include 
•^He in the Monte Carlo, we add scattering functions that 
account for the optical depth, mass difference and isotope 
shift in ^He. Photons are injected to the Monte Carlo in 
either the ^He or ^He lines, i.e. their initial frequency 
distribution is the sum of two Voigt profiles with nor- 
malizations in the ratio /4('^He). A photon is considered 
to have escaped the line when it either redshifts out or 
photoionizes H i. 

We track a combined ^Hc/'^He occupation fraction and 
escape probability for both isotopes in the level code. In 
principle, one should use the Monte Carlo to develop a 
2x2 matrix of probabilities for absorption of photons by 
'^He/'^He given that it was emitted in an incoherent pro- 
cess in ^He/'^He, and then develop individual rate equa- 
tions for "^He and ''He occupations in the level code. We 
expect that this would be a small correction because in- 
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FIG. 6: Comparing the effect of Thomson scattering, ^He 
scattering, and feedback in the forbidden and allowed lines. 
The uppermost curve is the difference between two models 
with and without Thomson scattering in the n^P° — l^S and 
n^D — l^S where neither has feedback. This retards recombi- 
nation because, on average, more photons are injected into the 
optically thick region of the line. (Note that once feedback of 
the radiative distortion is added, the effect of Thomson scat- 
tering is greatly reduced.) Thomson scattering in 2^P° — l^S 
decreases the escape rate at early times, but once H I opac- 
ity becomes significant, more photons are on average removed 
from the optically thick part of the line. '^He can be seen to 
accelerate recombination slightly by assisting photons out of 
optically thick regions at late times. We note, though, that 
the typical error induced by resampling the Monte Carlo is of 
order 10^'* (see Fig. I1U|I . To confirm the effect of ^He would 
take a significantly finer grid of probabilities, but the actual 
value is immaterial to He I recombination overall, if it is this 
small. 

eluding '^He in the first place has a negligible effect; how- 
ever this assumption should be revisited in future work. 

Fig. [S] shows the cumulative effect of Thomson scatter- 
ing of Sec. Inland '^He scattering. The variation in Xe from 
■^Hc scattering is negligible and of order 2x 10~^, which is 
of the same order as the error induced by resampling the 
escape probability grid - indeed, to confirm that there is 
any effect at all would require significantly more Monte 
Carlo information. 



IV. ADDITIONAL EFFECTS AND SOURCES 
OF ERROR 

A. Rare decays 

In Paper I we considered the effect of feedback and con- 
tinuous opacity in the P" -1^ S , 'n?P°-l^S, n^D-l^S 
transitions in "^He. The spontaneous transition rates for 
n^D - 1^5 m and the allowed [H, [13, [H, [3§, ^ lines 
are well-known, but there is considerable variation in 
the rates described in intercombination literature. Three 
accurate methods of finding the intercombination rate 
exist, namely the relativistic random phase approxima- 
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tion [11], H^lleraas-type [U methods and the Z- 
expansion . We used the rates from Ref . [1^ as these 
were available for the greatest number of states. 

We have not resolved the fine structure J-levels in 
ortho-He i in our multilevel atom code, assuming instead 
that each level is populated according to its statistical 
ratios. This assumption is valid because the only rates 
that depend on ,/ (instead of merely the quantum num- 
bers n, iS*, and L) are the intercombination line rates, 
and these are all slow compared to allowed transitions in 
ortho-He i that mix different values of J. For example, 
mixing among the i ^ levels by emission followed by 
absorption of a 1.08 /im photon is ~ 10'"' times faster than 
the intercombination decay "^P^ \^S. The fine struc- 
ture also does not affect the Voigt profile of the He i] 
23po — l^^ line because, on account of angular momen- 
tum conservation (which is exact) only the J = 1 fine 
structure level can participate in an electric dipole line 
connecting to l^S". In the He i] 2^P°-\^S literature, the 
quoted spontaneous decay rate [1^ is for J = 1, which 
accounts for only 3 of the 9 possible states in the 1?P° 
terms that arc not resolved in the level code. Thus, the 
spontaneous rate used here, ^23p°^iis is 1/3 of the typ- 
ical literature values. 

Fig. [7] considers the effect of variation in the 2''P° — 1^5* 
spontaneous rate, the contribution of intercombination 
transitions from n > 2 (which are taken from [2^ and 
scaled by rT^ for higher levels), and the vi^D — se- 
ries. These uncertainties are shown to be considerable 
for He I recombination, and outweigh many more subtle 
effects. (Also note that the magnetic-quadrupole transi- 
tion 2^P2° ~ has a spontaneous rate ^ 0.4 s~^ [i^ . 
and is negligible.) 

There is also a metastable triplet level 2^S\ with a for- 
bidden one-photon decay, but its rate is very small, A = 
1.27 X 10~^s"^ [11]. This line has not been included be- 
cause it does not contribute significantly to the formation 
of ground-state He i. A rough estimate of its importance 
relative to decays from the metastable singlet level can 
be obtained as follows. The ratio of level populations will 
be roughly Sg^-E^/'^bT, (^here AE/kn = 9240 K) since the 
excited levels are in equilibrium during He i recombina- 
tion to a very good approximation. The ratio of rates of 
formation of ground-state He i via decay from 2^5 versus 
2^S will then be - 3e'^-^/'=BT,(^^^,.pj^^/A^.^gj^^)^ which is 

< 2.2 X 10"'' at z > 1000, which is negligible. 

Another route to the ground state is the electric oc- 
tupole decay [He i] n^F"-l^S with n > 4. The rates for 
these have not been calculated to our knowledge, however 
generically one expects them to be suppressed relative to 
the electric quadrupole transitions by a factor of 0{a^). 
Given that the quadrupole decays produced a change in 
the ionization history of order |Axe|max ~ 3 x 10~^, and 
given that the octupole decays come from higher-energy 
levels than with correspondingly lower abundances, 
we expect a smaller effect from the octupole series (the 
effect need not be a factor of ~ less because the [He i] 
3^D-1^S line is optically thick during much of helium 
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FIG. 7: There is significant disagreement in the literature 
describing the intercombination rate because of the approxi- 
mations involved. Here we show that the influence of a ±50% 
difference in the 2'^P° — l^S spontaneous decay rate (here 
denoted A2) on Axe is roughly < 1 x 10"'^ relative to the 
rate used throughout of 57 s~^. The uppermost plot shows 
the case where the 2^P° - l^S rate is neglected entirely. 
Once the feedback of non-thermal distortions between lines 
is taken into account, the overall effect of the intercombina- 
tion lines becomes much less significant. (Note that all the 
escape probabilities used here include continuous opacity, de- 
veloped in Paper I). Also included are the cases where the 
entire n'^P° — l^S series is included, and the n^D— 1^ S series 
are included in the base model. The n > 2 intercombination 
rates are truly negligible, and the contribution from n} D — l^S 
is of order, or less than the uncertainty due to the 2^P° — l^S 
rate. 

recombination). In fact, the effect of including octupole 
transitions is probably much less than quadrupole transi- 
tions because they overlap: the 4^F°-1^S' and 4^D-1^S' 
lines are split by Avjv = 2.8 x 10^^, which is 2.5 Doppler 
widths at z = 2000, and the splitting is less for higher n. 

B. Collisional effects in He I 

Collisions of atoms and ions with electrons play a key 
role in determining level populations in the interstellar 
medium. In general, they are less prominent during cos- 
mological recombination because of the very high den- 
sity of photons. There are four major collisional pro- 
cesses that affect the recombination history: (Ibb) re- 
distribution of excited levels of H i. He i, or He ii; (Ibf) 
collisional ionization or its inverse, three-body recombi- 
nation; (2) excitation/dc-excitation of an atom from/to 
the ground state; and (3) charge transfer reactions. We 
argue that these considerations can be neglected for He i 
recombination. 

Here the discussion is limited to collisional processes 
in He i during He i recombination. Because collisional 
processes generally draw species into equilibrium, and H i 
is very nearly in equilibrium during He i recombination, 
collisional processes in H i can be neglected here. Moving 
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FIG. 8: The departure of several levels from equilibrium val- 
ues during He I recombination. The excited states begin to 
fall out of equilibrium for z < 2000. Significant among these 
are 2^P° and 2^5", departing from Saha values by of order 
10% and larger hy z = 1700. (Note that the departure his- 
tories of 2^P° and 2^5 are so similar that here we only show 
2^P°. Likewise with 2^5 and 2^P°, so only 2^P° is shown.) 
A significant concern is that coUisional rates could drive 2^P° 
and 2^S closer to equilibrium and modify the overall 1^ S for- 
mation rate. 



in to H I recombination, especially in its late phases, H i 
coUisional processes are extremely important as the levels 
fall out of equilibrium and / and I' sublevels can be easily 
shuffled by coUisional processes H^]. Late time effects 
in H I recombination deserve much attention, but are 
beyond the scope of this paper. 



1. Equilibrating 2^ S by rates connecting to n > 2 

During He i recombination, the highly-excited popu- 
lations in both He i and H i are very close to Saha- 
Boltzmann equilibrium. Fundamentally, rates for elec- 
tron collisions inducing bound-bound and bound-free 
transitions are small then because they depend on depar- 
ture from equilibrium. Fig. [8] compares the departure of 
several highly-excited states from their Saha occupations 
(here, without coUisional effects added, which would push 
the levels further into equilibrium). We note that 2^P° 
and 2^5* depart from their Saha occupations by tens of 
percent at z ~ 1700. CoUisional processes will tend to 
push these closer to equilibrium with the continuum, and 
this will modify the formation rate to l^S. We will see 
that the coUisional rates connecting 2^P° to excited lev- 
els are not sufficient to affect its departure from Saha. 

The electron-collision induced transition rates are 
given by the standard rate equations which respect de- 
tailed balance 



(23) 



where 



.9ne 



2TTmpkT„ 



3/2 



(24) 



for coUisional ionization and three-body recombination. 
For the coUisional bound-bound rates, we have 



coU,bb n'>n 



(25) 

Dn'^nne, (26) 



where 



^{E„-E^,)/{kTr^) 



g-a 



(27) 



coll^bf 



(Note that between the excitation and de-excitation we 
have reversed the subscript order Dn'^n to Dn,n', respec- 
tively.) The total rate n' ^ n to a level n then de- 
pends on the population of n', the departure of n from 
Saha equilibrium and the process's rate coefficient, here 
{av)ne, as in radiative bound-bound processes. We can 
thus compare the coUisional rates to the spontaneous rate 
Ai'-tn-Ps connecting the levels. For the allowed transi- 
tions 2^5* — n^P°, {av) is given by Ref. [i^: {av) ^ 
{5 X 10"^ 1 X 10-^ 4 X 10"^, 2 X 10-^}r4"^/^ cm^ s'^ for 
n = 2, 3, 4, 5 respectively. These are all between 9 and 12 
orders of magnitude smaller than the radiative Einstein 
coefficients connecting these levels (for Ue — 1080 cm^'^). 
The radiative rates connecting levels are not sufficiently 
suppressed by their escape probabilities for coUisional ef- 
fect to contribute (recall the allowed 2^P° — l^S series 
is optically thick, but transitions from and to n > 2 are 
thin to a good approximation). 

To assess whether coUisional processes through radia- 
tively forbidden lines can push 2^5* closer to equilibrium, 
we can compare the total coUisional rate into 2^5* to the 
total radiative bound-free rate supplying 2^S. If we con- 
sider decays from from n' to n in Eq. (j26[) , and define the 
departure from equilibrium as 6„ = a;„/n„ saha, the sum 
becomes 

■K{av)n'^nn^, (28) 
while for thermal radiation, the rate is 

Xn\bf = a::n,Saha(l - K) [^2^ S ■ (29) 

When 2^5" begins to fall out of equilibrium the typical 
coUisional rate coefficient {av)ne to all other levels [i^ 
is between 10"^ s^^ and 10^^ s^^. Combined with the 
factor exponentially suppressing transitions of increasing 
energy separation, it is clear that n = 2 dominates the 
sum. Yet, these rates are dwarfed by the radiative bound- 
free rates to the level, for which /32is ^ 10'^ around 
z ~ 1700. 
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2. De-excitatation of the n = 2 family to l^S 

It is possible to speed up recombination if collisions 
efficiently de-excite helium atoms to their ground state 
without producing further ionizing photons. In the case 
of He I there is the possibility that the collisional de- 
excitation of the metastable 2^5* level could compete with 
two-photon decay, which starts from the higher-energy 
2^5* level. Assuming approximate equilibrium between 
215' and 235*, n{2^ S)/n{2^ S) « Se^'^/^^^ « 25 at z ~ 
1600. The colhsion strengths of Ref. [13, SI give (av) - 
2 X 10~^T4~^^^ cm'^ s"^ where T4 is the temperature in 
units of 10"* K. The ratio of collisional de-excitation to 
two-photon de-excitation is then 



► lis 



(30) 



which is of order 10 ® throughout He i recombination 
and hence can be neglected. Considering the rest of the 
n = 2 family, (av) - {0.2,1,0.5} x crn^ s-^ 

for decays 2^P° ^ l^S, 2^S l^S and 2^P° ^ l^S, 
respectively. As in the case of 2^S, the relative occupa- 
tions in the n = 2 family are not sufficiently different to 
enhance overall collisional de-excitation rates for these 
processes to matter. Similar collision strengths for de- 
excitation to the ground state are obtained from higher 
levels of He i, but their occupation probabilities are much 
smaller, hence they will contribute even less to the for- 
mation of He I l^S. 

For He 11 recombination, the radiative rates are suf- 
ficient to maintain Saha equilibrium at all times to an 
accuracy of 0.2% in Xe (much less in Ci), and this cir- 
cumstance cannot be altered by including collisions. 



3. Charge transfer 

Charge transfer reactions will accelerate He i recombi- 
nation and proceed through 



He^ 



H ^ He + H+ + 7. 



(31) 



(The nonradiative reaction is much slower due to the 
large separation between potential energy surfaces for the 
initial and final states [4y|.) This is similar to continuous 
opacity, examined in Paper I, where a H i population can 
assist He i recombination. (Indeed, because continuous 
opacity acting within the line accelerates recombination 
as soon as a small H i population is present, H i popu- 
lations are not sufficient during He i recombination for 
charge transfer to occur.) The total charge transfer con- 
tribution to the rate is given by 



iHci = (1-2 X lO^^'^cm'^ s"') 



0.25 



300 K 



, (32) 



where = XHei — Xhi is the ionization potential dif- 
ference. Note that the second (detailed balance) term is 
only correct if the radiation and matter temperatures are 
the same, because it involves a photon; however this is 
the case throughout He i recombination. The rate coef- 
ficient is given by the fitting formula of Ref. [2§|. The 
charge transfer rate reaches a maximum at z ~ 1900 (if 
we turn it on in the interval z ^ 1700 to z ~ 2200), where 
|(ixHci/c?ln(l + z)| < 2 X 10"^. Thus it can be neglected. 



C. Peculiar velocities 

It has been suggested in the past that peculiar veloci- 
ties could be important for the recombination history of 
the universe. Taken at face value they could substan- 
tially affect the line profiles: typical baryonic peculiar 
velocities before matter-radiation decoupling are of or- 
der severalx 10~^c as a consequence of the primordial 
fiuctuation level, as compared with typical thermal ve- 
locities ^ (fcBTni/TOHc)^^^ ~ lO^^c. The similarity is (to 
our knowledge) a coincidence, but it implies that peculiar 
velocities could lead to significant broadening of the line 
profile. Seager et al. [43 neglected this possibility be- 
cause they treated line transfer using the Sobolev escape 
probability, which docs not depend on the line profile so 
long as it is narrow. In contrast, our calculation finds 
a very significant infiuence of the line profile on the re- 
combination history, so the issue of line broadening from 
peculiar velocities must be considered more carefully. We 
argue here that there is in fact not a major effect, at least 
in the standard cosmological model, because the peculiar 
velocities are coherent on large scales and hence the pe- 
culiar velocity gradients are small. Our argument may 
not apply to small-sc ale p eculiar velocities generated in 
more exotic scenarios j48l l49j , and it is even conceivable 
that recombination may provide a way of constraining 
these. Consideration of such exotic scenarios is beyond 
the scope of this paper. 

In some astrophysical applications, it is possible to 
treat peculiar velocities as providing an extra, non- 
thermal increase in the line width. This "microturbu- 
lence" limit applies when the mean free path is large 
compared to the coherence scale of the velocity field, and 
in this case it is only the distribution of velocities that 
matters (or only the rms velocity, if the velocity is Gaus- 
sian). The opposite limit is where the peculiar velocities 
have a long wavelength relative to the distance a photon 
can travel over the relevant timescale. In this case, the 
velocity gradient tensor diVj is constant over the region 
sampled by a photon, and it is instead the velocity gra- 
dient (in combination with other scales in the problem) 
that matters. One can see the difference between these 
two cases by considering what happens to a photon's fre- 
quency in the baryon reference frame: in the microtur- 
bulence case the frequency oscillates wildly around some 
mean value as the photon passes through eddies of dif- 
ferent baryon velocity, whereas in the uniform gradient 
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case the photon's frequency decreases hnearly, with the 
peculiar velocity gradient providing a correction to the 
Hubble constant (and providing it with direction depen- 
dence since diVj will in general have a spin-2 component). 
Since most of the velocity power during recombination 
is on scales of several Mpc comoving (the Silk damping 
length; see below), whereas the Thomson mean free path 
(anetTx)"^ is hundreds of kpc comoving (and even less 
for photons near resonance lines), microturbulence does 
not apply here and a more sophisticated analysis is nec- 
essary. We investigate the nature of the baryon velocity 
field in Sec. IIV C 1[ and then consider how far a photon 
can travel during the time it spends in a resonance line 
in Sec. IIV C 21 Most of the numbers describing the ve- 
locity field and photon transport here will be quoted at 
z = 2000, however the qualitative picture is the same at 
other redshifts during helium recombination. 



1. The baryon velocity field 

We compute the baryon velocity power spectrum P-^ (k) 
using the Boltzmann code COSMICS [HO] in Newtonian 
gauge. We assumed a primordial power spectrum with 
= 2.4 X IQ-^ at fc = 0.002 and = 0.95, consistent 
with the recent WMAP results [111 . We then obtain the 
velocity structure function at comoving separation r. 



(a) Baryon velocity power spectrum at z=2000 



52(r)^([v(0)-v(r)]2)= / 2[l-jo(A:r)] 



,k^Py{k) dk 



27r2 k ' 
(33) 

where jo is the spherical Bcssel function. The result is 
shown at z = 2000 in Fig. O Note the importance of the 
Silk damping scale k^^ = 2.7 Mpc (as computed using 
the formula in Zaldarriaga & Harari [5l|): perturbations 
at A; 3> fcD arc damped out, with the result that at sep- 
arations r <C fcj)^ the relative velocities of packets of 
baryon fluid are dominated by expansion and shearing 
due to long- wavelength perturbations {k^^ ^ r). In this 
limit we find 



where 



92 

^rms 



fc3p„(fc) dk 
27r2 k' 



(34) 



(35) 



is the rms peculiar expansion of the baryon fluid. At 
z ^ 2000 this is Ormsa = 2.9 kms^^ Mpc"^ (In prin- 
ciple the small-scale structure function depends on the 
shear as well as the expansion, however for scalar per- 
turbations they are not independent.) On scales small 
compared to k^^ the baryon velocity fleld essentially has 
a uniform gradient diVj (though its expansion rate may 
be perturbed and its shear nonzero); on larger scales the 
Universe contains many patches of different baryon ve- 
locity. 




k (Mpc" 

(b) Baryon velocity structure function at z=2000 




r (Mpc) 

FIG. 9: The peculiar velocity (a) power spectrum and (b) 
structure function of the baryons at z — 2000. The vertical 
lines show the Silk damping scale (/cd = 0.37 Mpc~^, or k^^ = 
2.7 Mpc in the lower panel) and the Hubble scale aH/c = 
0.007 Mpc"^ or c/aH = 143 Mpc. 



The velocity field at other redshifts of interest (1700 < 
z < 3000) is qualitatively similar, the major differences 
being that the Silk scale is shorter at higher redshift, 
and the acoustic scale (which sets the oscillations in the 
power spectrum) is also shorter at high z. The rms pe- 
culiar expansion declines from ^rmsO = 5.4kms~^ Mpc~^ 
at z = 3000 to 6l„„sa = 2.2 kms'^ Mpc'^ at z = 1700. 
It only undergoes qualitative changes after kinetic de- 
coupling (z ~ 1100) after which the baryons fall into the 
dark matter potential wells. This later era is not relevant 
to this paper, although it could affect the late stages of 
hydrogen recombination. 

Before we continue, it is worth noting that if Silk 
damping were somehow turned off, then the velocity fluc- 
tuations of order lOkms"^ seen in Fig. [HKa) would con- 
tinue out to arbitrarily large k. These small-scale fluc- 
tuations would be in the microturbulence regime, and 
because they are larger than the helium thermal veloc- 
ity, they would significantly broaden the He i lines. In 
Appendix [B] we describe cosmological perturbation the- 
ory relevant for radiative transport in helium in more 
detail, considering coupling of helium to the other fluids, 
damping; and plasma oscillations. 
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2. Photon transport 

We now ask how far a photon can diffuse during the 
time it spends in a resonance hne. By combining this 
information with the results from Sec. IIV C 1[ we will be 
able to determine the typical velocity differences within 
this region, and what effect (if any) they will have on re- 
combination. We consider the two lines most relevant for 
helium recombination - the allowed He I 2^P° — l^S line, 
in which a photon typically undergoes many scatterings 
during passage through the line, and the scmiforbidden 
He i] 2^P° — l^S line, which has optical depth of at most 
a few. 

For the allowed He i 2^P° — 1^5 line, we wish to 
know how far the photon can travel as it redshifts from 
v+ = fiine + Afiino to v_ = i^une - Az^iino. (Recall that 
Atiine is the detuning beyond which the integrated op- 
tical depth to incoherent processes is unity. We only in- 
cluded incoherent processes in this definition since coher- 
ent scattering does not change the ionization/electronic 
state of the atom.) Since we have /inc <SC 1, the co- 
herent scattering optical depth is very large even in the 
incoherent-thin part of the damping wings. Therefore 
the photon moves by diffusion, and we should calculate 
its rms diffusion distance. 



idiff — - 



Cimfp dt 



1/2 



dv 



Huiinedr/dv Offline 

-1 1/2 
dv ^ 



1/2 



aHvii, 



(36) 



(Here L^fp is the physical mean free path. The denom- 
inator in the last integral is dr/dv] note that this is 
an upper limit if Thomson opacity is also important.) 
The integral is dominated by the damping wings where 
(/)(z^) = riinc/47r^Az^^, so we have 



2^/2 7r 



cAf, 



3/2 



^diff 



V3 affl/linoTs^^rJ/^ 



„3/2 



(37) 



where in the last line we have substituted Ai'iino = 
riinc'''s/inc/47r^ (scc Papcr I). At z = 2000, we have 
Ts = 2.0 X 10^ /inc = 2.8 X 10-3, giving Ldiff = 3.1 kpc. 

For the scmiforbidden line, the damping wings are neg- 
ligible and the profile is Gaussian with width Ai^d- For 
most of the recombination history we cannot define an 
optically thick part of the line in analogy to Az^iino, be- 
cause the optical depth is ts < 2 and hence either half 
of the line is optically thin. We can however calculate 
the comoving straight-line distance a photon must travel 
before it redshifts by Ai^d, which (since the optical depth 



never exceeds a few) should serve as a guide to the typ- 
ical distance a photon can travel while interacting with 
the scmiforbidden line. This distance is 



aHviine 



(38) 



or 2.2 kpc at z = 2000. This is an upper limit on the net 
distance traveled by the photon because the latter could 
be reduced by scattering (although for optical depths of a 
few this reduction is likely to be modest). Similar results 
hold at other redshifts: at 1700 < z < 3000, we find that 
Lsi has a peak value of 2.3 kpc, and Ldiff has a peak of 
4.0 kpc. 

It is easily seen that the distances Ldiff and Lgi are 
much less than the Silk length. This is not surprising: 
the Silk length is (roughly) the distance a photon can 
diffuse during a Hubble time, whereas we are asking how 
far a photon can go during the much shorter time re- 
quired to redshift through a line. These distances corre- 
spond to scales over which the rms velocity difference is 
S'2^^(Ldiff) = 9ms^^ and ^^^(Lsi) = 6ms-^. In com- 
parison, the rms thermal velocity per coordinate axis is 
(^sTm/TOHe)"'^^^ = 3.4 km s"^. Thus over the region that 
a photon explores during its transport through a line, the 
peculiar velocity differences are very small compared to 
the thermal velocities of the helium atoms. This is true 
even though the baryons have a bulk (>Mpc scale) flow 
at an rms velocity of 18kms^^. 

Even though the relative peculiar velocities are small 
compared to the thermal velocity, it is still possible that 
the velocity gradient can lead to a net effect on the pho- 
ton's frequency after many scatterings. In order to assess 
this possibility, we measure the photon's frequency in the 
baryon frame rather than in the Newtonian frame, as this 
leaves all scattering terms in the radiative transfer equa- 
tions unchanged. It also eliminates the explicit depen- 
dence of position on the problem, so long as the velocity 
gradient is uniform over the region explored by the pho- 
tons (which is true since Lsh^diff k^^). In this case, 
the photon is moving through a medium whose expansion 
rate (i.e. velocity gradient) is not exactly given by the 
Hubble rate H, but rather by a correction H+9b/3, where 
Ob is the expansion perturbation. The velocity gradient 
may also have a shear, which for scalar perturbations is 
generally of the same order as 9i,. The Hubble expansion 
term i'l Hubble = —Hv in the radiative transfer equation 
is replaced: 



(39) 



where is the baryon shear tensor and is the direc- 
tion of propagation. The peculiar expansion has vari- 
Since we have for each scalar Fourier 
it follows that 



ance 
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rms 



5„-/3)06(k), 



mode cr{,ij(k) = {kikj 
{<Jb,ij<^b,ij) = (2/3)0?^,. 

The effect of the 9b term in Eq. ([M)) is equivalent to a 
local modification of the Hubble rate, with rms fractional 
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change O^^s/'^H = 5 x 10~^. Its effect can therefore be 
modeled by changing the Hubble rate and asking what 
effect this has on the recombination rate xhcI- The ef- 
fect of the Hubble rate is to change the escape probabil- 
ities, whose rms change will be at most 5 x 10"'' (escape 
probabilities scale as Pg oc i7 for optically thick lines 
with no continuum opacity, and the scaling is shallower 
if there is continuum opacity or if rg is of order unity or 
less). Thus the typical correction to the recombination 
rate (5iHci/iHoi will be of order 5 x 10~^. This correction 
will be positive in some regions and negative in others, 
because {Ob) = 0. Because ab has spin 2, the change in 
the recombination rate ^iuei can only exist at order ct^, 
and will probably be much smaller than the effect of the 
Ob term. 



tons through photo-detachment/photoionization. (neg- 
ative hydrogen in an important part of the opacity in 
stellar atmospheres [H^) Taking a cross section for neg- 
ative hydrogen of « 10~^^ cm^ [s^ and assuming Saha 
populations, the differential opacity is « few x 10^^ times 
lower than H I (near the end of He i when the negative 
hydrogen depth is greatest), being largely suppressed by 
the Saha occupation for the binding energy of 0.755 cV. 
While the ionization threshold of lithium is 5.39 eV [soj . 
its fractional abundance is at the level of 10~^ — 10~^° 
[60t . Using the cross sections to He i 2^P° — radia- 
tion from [6l|, we find that the differential optical depth 
to lithium, is a further three orders of magnitude smaller 
than negative hydrogen. Both processes are truly negli- 
gible. 



D. Additional considerations 

There are several issues that we have neglected 
throughout. The first is the modification to the vac- 
uum dispersion relation from multiple resonant scatter- 
ing, such as through 2^P° — l^S :52;]. (This is associated 
with a modification to the refractive index, or alternately 
the complex dielectric constant and group velocity, in the 
neighborhood of the line.) Throughout, we have assumed 
Lu = ck and that the radiation fields arc completely de- 
scribed by Af{uj) in calculation of the matrix elements. 
Now both of these assumptions break down, and the issue 
arises that the spontaneous and stimulated decay rates 
are modified by the mild dielectric behavior of the gas 
near the resonance frequency. These rates are the subject 
of active research [HI, [13, HM, [H, [111 . Generally, the mod- 
ifications are directly related to the real part of the refrac- 
tive index, the group velocity, and the imagin ary part of 
the dielectric constant (from absorption) [5J, [50| . These 
have maxima U{n) - 1 w 6 x 10"", 9(e/eo) « 2 x 10'^^, 
1 — Vg K. 7 X 10~^^ around the He i 21.2 eV resonance. 
A further treatment would require a significant overhaul 
of the methods used throughout, and because of its sub- 
tlety and complexity (how does radiative transport work 
when the radiation and matter statistics arc correlated; 
what is the appropriate density of states; how can this 
be implemented economically?) would vastly complicate 
the argument. We believe that because the pertinent 
quantities that describe deviations from vacuum propa- 
gation are so small that this is a negligible addition to 
the model's physics. There are lower-frequency plasma 
effects associated with the plasma frequency (typically 
several hundred kHz here), which leads to negligible dis- 
persion for any of the resonances involved, especially be- 
cause the relevant quantities depend on uj^juj^ for some 
transition frequency loq. 

In paper I, we emphasized the corrections in transport 
through the He i line due to a small population of H i. 
The negative hydrogen ion, neutral lithium [H^l, HeH^, 
H2 and H+ are also present, albeit with very low occu- 
pations, and can act as a sink for He i resonance pho- 



An upper limit to the effect of HJ and HeH+ can be 
obtained by assuming that their abundance is in ther- 
mal equilibrium according to the reactions H(H"'',7)H^ 
and He(H+,7)HcH+ (there are other pathways for for- 
mation of these species, but they are suppressed here 
by the abundance of the precursors and the tempera- 
ture), and that their cross sections to UV photons are 
given by the unitarity bound for electric dipole transi- 
tions, a < SA^/Stt. We note here that the equilibrium 
abundance of these ions is rapidly increasing as the uni- 
verse expands and cools, so the true abundance may 
be much less than this if the reaction rates are slow 
(see e.g. [S^)- Also photoexcitation and photodissoci- 
ation cross sections are usually several orders of magni- 
tude less than the unitarity bound except at the cen- 
ters of resonance lines, so the analysis here is proba- 
bly very conservative (i.e. overestimates the effect by 
a large amount). The Hj levels used here are the same 
as those computed in Ref. [63 | (with a binding energy 



taken to be 2.65 cV [62 
tion is taken from Ref. 




and the HeH"*" partition func- 
631 1 (with a binding energy taken 



to be 1.84 eV [64|.) For H2 , the equlibrium analysis gives 
x[Hj] 4 X 10^21 at ^ = 1700 (the end of He i recombi- 
nation), for which the unitarity bound implies a depth of 
nac/H < 5 X 10~^ at 21.2 eV; this number is much less at 
higher redshift. For HeH+, we find a;[HcH+] - 4 x 10-^5 
at z = 1700, hence nac/H < 5 x 10"^". Again this num- 
ber is much less at higher redshift. For H2, using the 
partition function from [ssj (assuming equilibrium pop- 
ulations) we get a;[H2] ~ 2 x 10"^'^ &t z = 1700, giving 
n<7c/H < 2.8 X 10~®. Thus opacity due to these species 
will have no significant effect on He i recombination. 



Also neglected is the possibility of a velocity-dependent 
occupation fraction for the atomic levels. We have 
treated the atom's absorption profile conditioned on the 
struck atom's velocity (which is taken to be thermal) 
in the Monte Carlo developed in Paper I. We have ne- 
glected these effects in incoherent scattering by assuming 
the photon is re-emitted across a pure Voigt profile at line 
center. These corrections are expected to be small, but 
warrant consideration. 
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E. Convergence of the methods 

Here we check the convergence and accuracy of the 
numerical methods employed throughout this series of 
papers. The most significant among these for the overall 
rate are the level of refinement of the probability grids 
estimated in the Monte Carlo and their accuracy (which 
can be assessed easily by resampling) and the accuracy 
of the atomic level code numerical solution. 

Here we consider three refinement and resampling 
cases: 1) doubling the number of MC photons in the 
sample in the 11 x 21 (redshift by xhci) grid of probabil- 
ities, 2) resampling the 11 x 21 grid, and, 3) refining to 
a 21 X 41 grid. These arc shown in in Fig. [101 Wc also 
confirm that the feedback iterations have converged (in 
the sense of giving negligible differences between subse- 
quent iterations) in Fig. [TlJ We also considered the level 
code step size. By doubling the number of steps through 
the recombination history, wc change the recombination 
history by |Aa;e| < 3 x 10"^. 

In this series of papers, we have ignored levels with 
n. > fT-max- As the matter and radiation temperatures 
drop, n„iax must increase to account for levels that have 
possibly fallen out of equilibrium. If the net recombina- 
tion (capture) rate to these high states becomes signif- 
icant, then truncation at n„iax will generally contribute 
less to the formation rate of the neutral species. Here 
we simply (roughly) halve rimax to rimax = 45 (from 
?^max = 100), to find that the change in Xe{z) is neg- 
ligible and of order \Axe\ < 4 x 10~^. The contribution 
to the formation rate of the ground state from the de- 
cay of these highly excited levels is greatly suppressed 
by the feedback of the spectral distortions they generate. 
Indeed, given the high optical depth in the n^P° — 1^5* 
lines, a model neglecting feedback over-estimates the con- 
tribution of these highly-excited states to Xe by slightly 
over an order of magnitude. 
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FIG. 10: Comparing several numerical convergence issues in 
the Monte Carlo-estimated 11 x 21 (redshift by a;Hoi) grid 
of escape probabilities. Both doubling number of photons in 
the sample and resampling the Monte Carlo give corrections 
of order < 2 x 10~*. Grid refinement is a more significant 
systematic, roughly < 4 x 10"'*, and indicates that log-log 
interpolation on the coarser grid over-predicts the escape ve- 
locity. Halving the step size in the level code results in error 
lAsel < 3 X 10^^ 
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V. EFFECT ON CMB ANISOTROPY 

Our principal motivation for developing a detailed un- 
derstanding of He I recombination in this series of papers 
has been to improve the accuracy of CMB models. This 
section and the following section are devoted to the mag- 
nitude of improvements to the CMB power spectrum Ce, 
and to the residual uncertainty from He i recombina- 
tion. Wc will argue here that the residual uncertainty 
is small compared to the dispersion among some of the 
recently discussed recombination histories. The residual 
uncertainty from He i recombination in the CMB spec- 
tra still exceeds cosmic variance (because of the large 
number of modes available) and is comparable to the un- 
certainty from the precision of modern linear Boltzmann 
codes [661 . This uncertainty is effectively confined to a 
one-parameter family of {Ci}s parameterized by the Silk 
damping scale. 

The helium recombination history affects CMB 



FIG. 11: Convergence of the iterations to include feedback of 
non-thermal distortions between lines. These are descending 
from the diflference between no feedback and one interation, 
between the first iteration and the second, and so on. Note 
that by the 4th iteration, the efi^ect is roughly Axe < 10~*, so 
by going a fifth iteration, any systematic effect is negligible. 
Note that the integration tolerance taken in the level code is 
1 X 10~^ 



anisotropics by changing the electron abundance in the 
redshift range 1500 < z < 2700. In general the change in 
CMB power spectrum can be written as 



ACT 

(jTT 



j FT{z)AXe{z) - 



dz 



(40) 



where the function Ft is written in terms of the func- 
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tional derivative, 



Fj{z) = 



1 + z 

cr 5xe{zy 



(41) 



That is, an increase of Axg = 0.01 for N e-folds of 
expansion results in an Fj N% change in the temper- 
ature power spectrum, so Ff {z) can be considered as 
a weighting function that indicates how much Cj'^ de- 
pends on a given point in the recombination history. For 
polarization-sensitive experiments one may also define 
an analogous function Ff{z). Since the temperature- 
polarization cross-spectrum Cf^ crosses zero twice every 
acoustic oscillation, we choose to define 



i^/(z)==(l + z) 



TE 



SXe{z) ' 



(42) 



where pj^ = Cj^ / y Cf^Cf^, so that the temperature- 
polarization correlation coefficient is modified by 



TE 



Fx{z)Axe{z) 



dz 
l + z' 



(43) 



To reach sampling variance accuracy over some range 
Ai, each of ACf^/Cf^, ACf^/Cf^, and pf^ must 
be calculated to an accuracy of order 1 / V(A£, which in 

the damping tail is of order 3 x 10^** /^j^^^^. Since he- 
lium recombination lasts for roughly 0.5 e-folds, it fol- 
lows that the accuracy requirement on Xe{z) is roughly 
6 X lO^'^F^^/gj^.y^'^. In practice, foregrounds (mostly ex- 
tragalactic) and beam uncertainty will limit the accuracy 
of small-scale CMB experiments and "Fisher matrix" ac- 
curacy on the high multipoles may not be reached, nev- 
ertheless this docs provide a target for the theoretical 
calculations. 

The functions F^'^'^{z) are displayed in Fig. [T^ We 
have measured these by numerical differentiation of the 
power spectra computed using the CMBFAST code [13 ■ At 
each redshift, a Gaussian of width cr^ = 100 in redshift 
and normalization J Axe dz = 50 was artificially added 
to the recombination history, and the change in Ci was 
used to estimate fJ''^''^ {z). (This is equivalent to tak- 
ing the functional derivative for sharp test functions and 
convolving by a Gaussian, and ensures the stability of the 
CMB model calculation.) Fig. [T21 shows several major ef- 
fects. One is the trend that F'[{z) rises from near zero at 
low multipoles to some positive value of order 1-3 for high 
multipoles, and is most significant for corrections to Xe{z) 
during later periods of He i recombination. This is sim- 
ply the result of a change in the Silk damping scale: the 
low multipoles are determined by low-fc Fourier modes in 
cosmological perturbation theory, for which the baryon- 
photon fluid is nearly ideal at these redshifts. The high 
multipoles are of course Silk-damped, and the sign of F 
is positive because increasing the electron abundance re- 
duces the photon mean free path and hence reduces the 
diffusion length, thus there is less suppression of power 



and Ce goes up. The variation in Xe{z) becomes more 
important at late times in He i recombination because 
of the increasing overlap with the visibility function. In- 
deed, this begins to give an indication of how important 
a solid understanding of H i recombination (which begins 
z - 1700) will be. 

For polarization, we see that Ff{z) w F'[{z), because 
the suppression of acoustic oscillations by Silk damping 
has roughly the same fractional effect on CMB tempera- 
ture and on polarization. (While it is true that some po- 
larization is generated during helium recombination due 
to the finite photon mean free path and consequent pro- 
duction of a quadrupole moment, the high Thomson opti- 
cal depth at that epoch guarantees that this polarization 
is erased before it reaches the observer.) This is also why 
the dimcnsionlcss correlation coefficient depends very 
weakly on recombination: |i^/'^(z)| <^ \fJ'^{z)\. 



VI. SUMMARY OF RESULTS FOR HELIUM, 
DISCUSSION 



In the series of Papers I-III, we have addressed a multi- 
tude of effects in He i recombination, some new, and some 
that have appeared recently in literature [l^, [lO, [HI, H^l ■ 
Here, we described the effect of Thomson and ^He scat- 
tering, peculiar velocities in the recombination plasma, 
and coUisional and rare processes. Table U summarizes 
the processes described in Papers I-III and the magni- 
tudes of their effects, and Fig. [13] shows the cumula- 
tive difference in the CMB temperature and polarization 
anisotropics due to the largest corrections described. We 
want to emphasize only the modification to the He i re- 
combination history from the effects described in Papers 
I-III, so we only compare with a reference model where 
those effects are absent. A summary of the final He i 
recombination history developed here can be found in 
Fig. [Ml We have not described He ii recombination in 
detail because it is almost a pure Saha rccombinaton, 
as shown in Fig. [15] and is insignificant for the CMB 
anisotropy. 

Thomson and '^Hc scattering contribute to Xe only at 
the level of a few x 10~^ and can be neglected (for com- 
parison, the uncertainty associated with the interpolation 
of the Pesc derived in Monte Carlo used here is only of 
order 10~^). The largest (known) systematic error pre- 
sented here is from plausible uncertainties [1^, [13, [ll] 
in the 2^P° - 1^5 rate and is of order |Aa;e| ~ 10"^. 
Making no assumptions about the correlations between 
systematic errors, a conservative estimate of the overall 
systematic in Xe{z) from known effects is 3 x 10^"^. 

The remaining uncertainty in the C^s due to helium 
recombination can be estimated by multiplying the max- 
imum error in the electron abundance, jAxelmax = 0.003, 
by / \FJ{z)\ dln(l + z) (see Eq. [41]). Over the redshift 
range of helium recombination, 1700 < z < 2800, we find 
/ \Ff{z)\ d\n{l + z) = 0.6 at £ = 3500, implying that 
Cj'^ has a fractional uncertainty of 0.002 (this number 
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TABLE I: Summary of the magnitude of effects described in Papers I, II and III. From the top category to the bottom, we 
distinguish the magnitude of an effect (the systematic error of not including it) from the uncertainty in an effect, and the 
uncertainty in the implementation. Unsigned upper bounds on \Axe\ are indicated by ±. Note that these are only meant to 
give an order-of-magnitude bound on the effect, where more detail is available in the section cited. The effect of forbidden 
processes includes H I opacity in their transport and feedback. The direction shown is either + (increases a;e) or — (decreases 
Xe), and a ~ indicates that only an order-of-magnitude was necessary to find that the effect was negligible. 



Effect 


Dir. 


1 max 




Sec. 


z 


Systematic corrections due to effects: 


opacity within lines (2^P° - 1^5) 




2.5 X 10" 


2 


I, HIE 


~ 1800 


feedback between 2^P° - l^S and 2^P° - 1^5" 


+ 


1.5 X 10" 


2 


I, nc 


1800 - 2600 


continuum opacity modification to feedback 




5 X 10" 


3 


I, IIIB 


~ 1800 


2^P° - l^S inclusion 




3 X 10" 


3 


in,|ivA| 


~ 1900 


n^P° -l^S,n>?, inclusion 


± 


~ 4 X 10" 


5 


in,|ivA| 


~ 2000 


tt}D — 1^5 inclusion 


± 


~ 3 X 10" 


4 


ni,|ivA| 


~ 1900 


opacity in P° - l^S and n^P° - l^S (for n > 3), n^D - 1^5 




5 X 10" 


4 


I, HIE 


~ 1900 


coherent scattering in 2^P° — l^S 




2 X 10" 


4 


I, HIE 


~ 2000 


distortion, thermal stimulated two-photon effects 


+ 


4 X 10" 


5 


II, II 


2000 - 3000 


electron scattering 


± 


3 X 10" 


4 


HI, El 


1800 - 2800 


Uncertainty in the effect's magnitude: 


finite linewidth in He I 


± 


~ 4 X 10" 


4 


II, vc 


1800 - 3000 


non-resonant two-photon effects from n > 2 


± 


~ 5 X 10" 


4 


II, IV 


~ 2000 


±50% 2^P° -l^S spontaneous rate 


± 


~ 10" 


3 


m.livAl 


~ 1900 


Uncertainty from the numerical implementation: 


modified escape probability grid refinement 


+ 


~ 5 X 10" 


4 


III|IVE| 


~ 1900 


Monte Carlo resampling 


± 


~ 2 X 10" 


4 


ra,|ivE| 


~ 1900 


Monte Carlo sample size doubling 


± 


~ 10" 


4 


in,|ivE| 


~ 1900 


level code: half step size 


± 


~ 3 X 10" 


5 


m,|iVE| 


~ 1800 


convergence in the modified P° — l^S series modified P^sc 




~ 3 X 10" 


4 


I, HI E 


~ 1900 


Timax = 45 relative to rimax — 100 


+ 


4 X 10" 


5 


III, IV E 


~ 2300 



TABLE II: Numerical summary of He I recombination in a model with continuous opacity and feedback between levels. We 
provide this for comparison with external recombination codes. The cosmological parameters used throughout are, Os = 
0.04592, Q.M = 0.27, Q.r = 8.23 x 10"^ Tcmb = 2.728, h = 0.71, /hc = 0.079 and zero curvature. 



Z Xe 


Z Xe 


Z Xq 




Z Xe 


3087.5 1.0790 
2996.2 1.0789 
2936.9 1.0788 
2878.7 1.0787 

2821.7 1.0783 

2765.8 1.0778 
2711.0 1.0769 


2657.3 1.0757 
2604.7 1.0743 
2553.1 1.0728 
2502.5 1.0711 
2452.9 1.0694 
2404.3 1.0677 
2356.7 1.0658 


2310.0 1.0639 

2264.2 1.0617 
2219.4 1.0594 
2175.4 1.0569 

2132.3 1.0539 

2090.1 1.0505 
2048.7 1.0464 


2008.1 1.0413 
1968.3 1.0349 
1929.3 1.0274 
1891.1 1.0196 
1853.6 1.0121 
1816.9 1.0057 
1780.9 1.0015 


1745.6 0.9999 

1711.0 0.9994 

1677.1 0.9988 



is lower for smaller £) . The same calculation for Cf^ also 
gives 0.002, and for pj^ we get < 3 x lO""* at £ < 3500. 
The latter number is roughly equal to the cosmic variance 
limit ^ ^max- Therefore we believe that our present cal- 
culation of helium recombination is sufRcient to predict, 
at £ < 3500, the temperature and i?-mode polarization 
power spectra to 0.2% accuracy and the correlation co- 
efficient pj^ to cosmic variance accuracy. Note that this 
uncertainty does not include the effects of hydrogen re- 
combination. Also, while we have made every effort to 



include all effects, we acknowledge that there could be 
additional processes that escaped our imagination. We 
encourage further work to ensure the completeness of the 
understanding of He i recombination. 

We conclude that a modification yielding more accu- 
rate He I recombination histories will give a negligible 
change in the WMAP results. However, the considera- 
tions presented here are significant for the next genera- 
tion of small-scale anisotropy experiments, which seek to 
measure the large multipoles at the percent level. The 
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FIG. 12: The functions FJ'^''^ [z) that describe how recombi- 
nation influences the CMB power spectrum, where individual 
contours are for the redshifts specified. 



changes shown in Fig. [13] would affect a cosmic variance 
limited experiment at the ^ la level at £max — 1500, 
and ^ 8(7 at £niax = 3000. This means that the effect is 
in princi ple at the ~ la level for the upcoming Planck 
satellite [631 (nearly so) for a high-resolution CMB 
experiment mapping 1-2% of the sky to imax = 3000, 
although it is possible that systematic uncertainties such 
as beam modeling and point source removal may pre- 
vent one from reaching "Fisher matrix" accuracy. Of the 
corrections considered, we propose that feedback of non- 
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FIG. 13: The cumulative efi'ect on temperature and polariza- 
tion anisotropies from H I continuous opacity and feedback 
during He I recombination, as calculated by CMBfast. We 
note that here we only consider the difference between the 
reference helium model with and without these effects. Com- 
parison of the full H I and He I history to standard methods 
such as Recfast will be the subject of later work. 



thermal distortions between the allowed lines, continuum 
opacity from H i photoionization, and the inclusion of 
the 2^P° — l^S [201 rate become part of the standard 
recombination model. We also emphasize that one of the 
largest uncertainties in the He i recombination treatment 
presented here is the 2'^P° — 1^5* spontaneous rate. In 
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FIG. 14: The He I recombination history from our multi-level 
atom code (solid line), compared to the Saha equation (long- 
dashed) and the commonly used three-level code recfast by 
Seager et al. [6^ (short-dashed). Both our analysis and that 
of Seager et al. find that He I recombination is delayed due 
to the n = 2 bottleneck. However we find a slightly faster 
recombination than Seager et al. due primarily to our inclu- 
sion of the intercombination line He l] 2"^P° — 1^5 and the 
accelerating effect of H I opacity. The latter effect causes our 
He I recombination to finish at 2 ~ 1700, whereas in recfast 
one-third of the helium is still ionized at that time. 



electron abundance near the peak of the CMB visibility 
function. (It will also affect the temperature evolution, 
21 cm absorption, and chemistry during the "Dark Ages" 
between recombination and reionization.) The calcula- 
tion of H I recombination and its effects on the C^s to 
sub-percent accuracy will be far more difficult than the 
calculations for He i presented in this series of papers. 
One reason is that there are 12 H atoms for every He 
atom, so to achieve the same absolute accuracy in Xe re- 
quires a much better accuracy in the ionization state of 
H than of He. Another reason is that the optical depth 
in the H i Lya line can be of order 10^ (instead of 10^ 
for the analogous 21.2 eV line of He i), which means that 
slow processes can play a larger role, and that the finite 
lincwidth effects could potentially be much more signif- 
icant. (The high optical depth also makes the photon 
Monte Carlo developed here computationally demand- 
ing for hydrogen.) Finally, the very high-71 states of H i 
become important at lower temperatures, and both l- 
dependent occupations and coUisional processes in these 
states will have to be carefully taken into account [i^ . It 
is because of these complexities that we have not solved 
H I recombination here; rather we plan to follow up this 
scries of papers with a similar scries on hydrogen in the 
future. 



He II recombination 




7000 



FIG. 15: The He 11 recombination history from the level code 
developed here. He 11 recombination is essentially irrelevant 
for CMB physics, and because of its rates it varies from Saha 
evolution at the level of < 0.2%. (for example, at £ < 3000 
the absolute difference in Cj'^ for the full model relative to 
Saha is < 3 X 10"^) 
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APPENDIX A: THE ELECTRON SCATTERING 
KERNEL CHARACTERISTIC FOR DIPOLE 
ANGULAR DISTRIBUTIONS 



In this section, we solve explicitly for the characteris- 
tic function of the electron scattering distribution over 
frequency (P(Ai/) taken to be in the Thomson limit). 



Paper II we show that the corrections due to two-photon 
[201 and finite linewidth effects are almost entirely neg- 
ligible. A fast, modified recombination code based on 
these corrections will be the subject of later work. 

While He i recombination is important cosmologically, 
it should only be considered as a first step toward a 
much broader treatment of cosmological recombination 
as a whole. This includes H i recombination, which is 
especially important for the CMB since it affects the free 



7(fc) 



ik{Aiy) 



(Al) 



This is used in the Fourier domain solution to the trans- 
port equations for complete redistribution with electron 
scattering and continuum opacity which is negligible over 
the width of the line, described in Sec. Ill Dl and is appli- 
cable for the intercombination and quadrupole lines for 
z > 2000. Substituting the variance ( = 2cr|)[l — cos(x)] 
and inserting the dipole angular redistribution function. 
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this becomes 



w{k) 



' D Jo 



4af 



4 + 2 



Ri{k) i?2(fc) 



2i?3(fc) 



(A2) 



where we have absorbed the Gaussian integrals i?2j 
and i?3 which are explicitly, 



Ri{k) = 


16 


^l-[2a|,fc4 


R2{k) = 


4 


1 - (2cr^fc2 


Rm = 


2 


(^1 _ e-^k^-l 



(A3) 



Note that for an isotropic distribution, the result can be 
expressed very concisely as w{k) = (4cr|j)~^i?3(fc). Note 
that these require caution numerically for small k because 
of the near-cancellation of several terms. This can be 
remedied by using a high-order expansion in small k (the 
first order cancels and leaves divergent terms in fc 0). 
We also note that while the relativistic electron scattering 
kernel is in general asymmetric about Av = 0, in the 
Thomson scattering approximation considered here the 
kernel is symmetric, leading to all-real Tu{k). 



APPENDIX B: FLUID DESCRIPTION OF THE 
BARYONS 



to charge separation, (3) collisions, and (4) photoioniza- 
tion/recombination, which can switch particles between 
the He and He"^ constituents. In principle one should 
check energy exchange rates (e.g. whether the electrons 
and ions thermalize to a common temperature), however 
since the baryon thermal pressure plays no role in Silk 
damping these are not relevant here. 

Siegel & Fry [t^] have considered electric fields due 
to charge separation and find that they efficiently pre- 
vent the electron and ion densities from departing signif- 
icantly from each other. In particular, during recombi- 
nation (z -fewxlO^) they find [See [zl, Eq. (24)] 



(Bl) 



in their analysis the ions (i) consisted entirely of protons 
but the addition of some He ii or He ill should cause no 
qualitative changes. (It should result in the ion expansion 
6i being replaced by the charge-weighted expansion of all 
the positive ions, Qjdjnj /ug, where qj is ion charge 
and Hj/ue is a number density ratio.) This equation 
results from a balance of the radiation pressure on the 
electrons, which is the driving term in their separation 
from the ions, with the electrostatic force that seeks to 
eliminate bulk charges. Here 9 is the peculiar expansion, 
ctt is the Thomson cross section, and p-y/pb is the ratio of 
photon to baryon densities. Plugging in the cosmological 
parameters gives 



In this appendix, we investigate whether the treatment 
of the baryons as a single fluid is adequate for the inves- 
tigation of the peculiar velocity field. In particular, we 
would like to understand whether the single fluid treat- 
ment accurately describes Silk damping. To do this, we 
have to determine whether the collisional relaxation rates 
are fast compared to the oscillation timescale at the Silk 
damping length, kjj^. We will show in this appendix 
that the collisional relaxation rates are indeed fast, so 
it should be accurate to describe the baryons as a fluid. 
Silk damping then proceeds as in the usual picture at all 
scales (ko ^ 0.37 Mpc"^) < fc < (fc^ ~ 2 Mpc~^) at 
z ~ 2000 (with smaller scales evolving in this range at 
earlier times [tO] ) . Below the photon free-streaming scale 
(the mean free path of photons), only electromagnetic in- 
teractions can influence the baryon velocity moment. We 
will consider the possibility of charge separation on these 
scales and show that it is a negligible contribution to the 
velocity structure. 

There are five constituents of the baryonic plasma dur- 
ing the recombination epoch: the electrons, ions (H+, 
He+, He^+), and neutral species (He). (During helium 
recombination there is very little neutral H, and its ve- 
locity structure is not relevant because its only signifi- 
cant role is to provide continuum opacity.) These species 
acquire and exchange momentum through several pro- 
cess [7l| : (1) radiation pressure, which is only signifi- 
cant for the electrons, by mass, (2) electric fields due 



where 



Atg = 6 X 10 



-20 , 



2000 



(B2) 



(B3) 



During an acoustic oscillation 9^ and 9i, arc similar, but 
both on the order of lOkms^^ or less. Thus so long as 
Atq is much shorter than the oscillation period, AtqUJ <^ 
1, the electron and ion velocities will be similar. 

Next wc come to the collisional momentum exchange 
times. Of interest are the electron-ion collision rate, the 
collision rates between ions of various species, and the 
neutral-ion collision rate. The collisional momentum re- 
laxation rate for particle type 1 against particle type 2 
(where we will take 1, and 2 to be electron-proton, and 
He ii-proton) for a thermal distribution and small differ- 
ential drift velocity is given by 



3V27r 



?^2 



■ ln(A) 



(B4) 



where A is the Coulomb logarithm, p 
is the reduced mass, and 



TO1TO2/ {mi+m2) 
VkeT/fi^. For roughly 



z > 6000, the doubly- ionized helium population exceeds 
the singly-ionized population. The momentum relaxation 
rate for doubly-ionized helium against protons is four 
times larger than the He ii-proton rate. For simplicity. 
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we will only consider the He ii-proton rate in Fig. [161 for, 
if it is sufRcient to relax He ii to the velocity structure 
of the other baryons, the He ill-proton rate must also be 
sufficient. Note also that the relaxation rate for He ii 
on protons has no dependence on the number density of 
He II, while the relxation rate for protons on He ii does. 

We will consider only proton collisions with singly- 
ionized helium, even though at for z > 6000, the The 
rates involving neutral He must be considered at z < 
3500 when He i is present; they are usually slower be- 
cause they lack the long-range nature of the Coulomb 
force. The dominant rate is that of resonant charge ex- 
change with He"*": 



He + Hc+ Hc^ 



He. 



(B5) 



The charge exchange momentum transfer rate is well- 
approximated by [7J| 

t.Hc+,Hc = 4.4x 10-i3nHe+(2T,n)i/2 

X (11.6- 1.041ogio(2r,n))s"\ (B6) 

where Tm is in Kelvins and nHe+ is in cm~'^. We have 
not considered momentum exchange between He and 
Hc+ due to photoionization/recombination, or He-proton 
scattering; if we did then this would only strengthen 
the conclusion that the momentum exchange rate to the 
charged fluid components is fast. 

The acoustic oscillation frequency at the Silk damping 
length is given by 



cfcr 



WD 



v/3(l + i?) 



(B7) 



where the damping wavenumber fco is obtained as in 
Ref. [in. The damping length runs from lO^^Mpc 
comoving at z = 2 x 10^ to 3 Mpc at z = 1600. (Before 
z = 2 X 10^ the usual computation is not valid because 
electron-positron pairs increase the opacity.) 

We have plotted the momentum exchange rates and 
compared them to the plasma frequency Wp and the 
acoustic oscillation frequency at the Silk damping length 
WD in Fig. [m The plot runs from z = 2 x 10* un- 
til z = 1600 when helium recombination has completed, 
for practical purposes. In all cases the relevant coUi- 
sional rates arc many orders of magnitude faster than 
the acoustic oscillations. This means that for wavenum- 
bers > 10~''Mpc we expect that hydrodynamics is 
valid. Silk damping will occur, and the exponential sup- 
pression of velocity perturbations predicted by the usual 
treatment is correct. 

A more detailed treatment is required in order to un- 
derstand what happens at scales below 10~^ Mpc comov- 
ing. Our physical expectation is that acoustic oscillations 
at such scales would also be damped by photon diffusion, 
with the Silk damping slightly modified by inclusion of 



positrons, and by the analogous process of neutrino diffu- 
sion (since these scales are well within the horizon when 
neutrinos decouple). Even if this docs not happen, the 
le+20 
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FIG. 16: Several scales in the recombination plasma relevant 
for helium peculiar velocities. The fastest here is the plasma 
frequency, followed by the electron-proton (and then He ll- 
p) momentum transfer rate. Note that the plasma frequency 
is significantly faster than both the Thomson rate and the 
frequency of an acoustic oscillation at the Silk scale. This 
greatly suppresses the magnitude of the charge separation, see 
Eg. IBll In the lower plot, we focus on the region neutral he- 
lium evolution {z < 3500). The charge exchange momentum 
transfer rates between He and He+ are much larger than the 
frequency of baryon acoustic oscillations at the Silk scale dur- 
ing the period of neutral helium recombination. This brings 
neutral helium into a common fluid with the other charged 
baryons. 



photon mean free path at 



2 X 10» is 10-^"' Mpc co- 



moving, so for k < lO^'^Mpc^^ the usual Silk damping 
calculation applies; since these scales are smaller than the 
Silk damping length they will be exponentially damped. 
For higher k it is less obvious what happens, but since the 
neutral He mean free path is much larger than 10~^° Mpc 
comoving during He recombination there can be no struc- 
ture in the neutral He peculiar velocity field at smaller 
scales. 
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